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Abstract 

Quantum versions of random walks on the line and the cycle show a quadratic 
improvement over classical random walks in their spreading rates and mixing times 
respectively. Non-unitary quantum walks can provide a useful optimisation of these 
properties, producing a more uniform distribution on the line, and faster mixing 
times on the cycle. We investigate the interplay between quantum and random 
dynamics by comparing the resources required, and examining numerically how the 
level of quantum correlations varies during the walk. We show numerically that the 
optimal non-unitary quantum walk proceeds such that the quantum correlations 
are nearly all removed at the point of the final measurement. This requires only 
O(logT) random bits for a quantum walk of T steps. 
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1 Introduction 

Quantum computing has produced a range of algorithms showing improve- 
ments over classical algorithms, among the most celebrated are Grover's search 
of an unsorted database [l| and Shor's algorithm for factoring large numbers 
In the search for new quantum algorithms, quantum versions of Markov 
chains are a natural place to look, since classical Markov chains provide the 
basis for some of the best known algorithms for hard problems such as ap- 
proximating the permanent of a matrix and fcSAT (with k > 2) fl. 
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Simple quantum generalisations of classical random walks spread quadratically 
faster on the line p, 0| , and mix quadratically faster on the cycle [3] • These 
promising early results were soon followed by several algorithms based on 
quantum walks. Shenvi et al [8] proved a quantum walk can solve the unsorted 
database search problem quadratically faster, and Childs et al [ij proved an 
exponential speed up for crossing a particular type of graph. Ambainis 13 



gives an overview of quantum walk algorithms, and Kempe [TJJ provides an 
introductory review of quantum walks and their properties. 

Quantum walks on simple one-dimensional structures remain a fertile test- 
ing ground for further research. It was shown numerically by Kendon and 



Tregenna [12[ and recently proved by Richter [13j, [14| that the addition of ran- 
dom noise or measurements to the quantum walk dynamics can optimise the 
spreading and mixing properties for quantum walks on both the line and the 
cycle. A detailed survey of the effects of such decoherence in quantum walks 



can be found in [15]. In this paper we examine how the interplay between 
quantum evolution and random noise or measurements produces optimal com- 
putational properties. The paper is organised as follows. In §2 we describe the 
basic quantum walk model we used for this study, and in §3 we extend this 
model to include non-unitary operations. We then outline how to count both 
quantum and random computational resources in a comparable way in §4. 
Section 5 introduces the measure of quantum correlations we are using, the 
negativity. In §6 we present our results for the quantum walk on the line, and 
in §7 our results for the walk on the cycle. Discussion and conclusions are 
given in §8. 



2 Quantum walks on the line and the cycle 



This study considers quantum walks taking place in discrete time and space, 
using a quantum coin to control the choice of direction. This is because we 
measure how "quantum" the walk is by examining the correlations between 
the quantum coin and the quantum walker. A different method for assessing 
"quantumness" would be needed for continuous-time quantum walks (intro- 
duced for quantum algorithms by Farhi and Gutmann [l6|). 

A discrete-time (coined) quantum walk dynamics consists of a quantum "coin 
toss" operation C, followed by a shift operation S to move the quantum walker 
to a new position. These are repeated alternately for T steps of the quantum 
walk, and the final position of the quantum walker is measured. For quantum 
walks on the line and the cycle, we have just two choices of which way to 
step, so the quantum coin is two-sided. We write \x, c) for a quantum walker 
at position x with a coin in state c G {+1, —1}- For a walk on the line, x G Z 
and for a walk on a cycle of size N, we have x G TL^. A classical random walk 



2 



can only occupy one location at any given time, but a quantum walker can be 
in a superposition of different locations. The full state of the quantum walk 
|\I/(t)) can be written as a combination of terms in each basis state \x, c), 

!*(*)) = E^cWkc), (i) 



where ip x ,c{t) £ C. By convention, the normalisation is defined to be 

EltM*)l a = i. ( 2 ) 

x,c 

When the quantum walk is measured (in the basis just defined), the walker 
is found in a single location with a definite coin state, but we cannot predict 
with certainty which state this will be, it could be any of the states in the 
superposition described by eq. (1). The probability of finding the quantum 
walker at position y with the coin in state b is given by 

P(y,b,t) = \(y,b\^(t))\ 2 = \^vAt)\ 2 (3) 
since the basis states are orthogonal, (y, b\x, c) = S xy Sb c - 

The coin toss and shift operators can be defined in terms of their action on 
the basis states \x, c), 

S\x, c) = \x + c, c) (4) 
C\x,c) = (\x, -c) +c\x,c))/y/2 (5) 
One can add more general bias or phase into the coin toss operation, see, for 



example, |17j], but this does not greatly change the basic properties of the 
quantum walk on a line or cycle, so we will consider only the unbiased case 
in this paper. For a quantum walk starting at position x = with the coin 
in a superposition state (| — 1) + i\+l))/y/2, (where i = y/—T) we can write a 
quantum walk of T steps as 



|vI/(T)) = (SC) J (|0, -1) + i|0, +l))/v / 2 (6) 

The solution for |\I/(T)) may be obtained by various methods such as Fourier 
analysis [H| and path counting [(J, and has been studied extensively. The key 
result is that spreading on the line proceeds linearly with the number of time 
steps. We can use the standard deviation ctq(T) of the probability distribution 
to quantify the spreading rate. For a quantum walk on the line, 

*q(T) = E * 2 P(x, c, T)^(^l-±y T, (7) 

asymptotically in the limit of large T. In contrast, for the classical random 
walk, the standard deviation is &c(T) = VT. 
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On the cycle, we are interested in mixing times rather than spreading. Mixing 
times can be defined in a number of different ways, we choose the definition 
given in Q, 

M(e) = min{T | V t > T : \\P(x,t) - P u \\ tv < e} (8) 

where P u is the limiting distribution over the cycle, and the total variational 
distance (TVD) is defined as 

\\P(x,T) - P u || tv = \P(x,T) - Pu\- (9) 

x 

A classical random walk on the cycle mixes to within e of the uniform distri- 
bution in time proportional to iV 2 log(l/e), where e can be chosen arbitrarily 
small. 

Pure quantum walks, on the other hand, do not mix to a stationary distri- 
bution. Their deterministic dynamics ensures they continue to oscillate indef- 
initely. There are several ways to obtain mixing behaviour, first explored by 
Aharonov et al Q. By defining a time-averaged probability distribution for 
the quantum walk, 

P(x,c,T) = ±J2P(x,c,t) (10) 

1 t=0 



they proved that P(x, c, T) does converge to a stationary distribution on a 
cycle, and that on odd-sized cycles the stationary distribution is uniform. A 
mixing time can be defined for P(x, c, T), 

M(e) =min{T | V t > T : \\P(x, t) - P u \ | tv < e} . (11) 

and Aharonov et al [7] proved that, for odd-sized cycles, M(e) is bounded 
above by 0(e~ 3 iV log N), almost quadratically faster (in N) than a classical 



random walk. Kendon and Tregenna 12|, observed numerically that M(e) 



0(N/e), this has been recently confirmed analytically by Richter [1 

Notice that we pay a price for our time-averaging: the scaling with the precision 
e is now linear instead of logarithmic. Aharonov et al [7] provide a fix for this 
in the form of a "warm start". The quantum walk is run several times, each 
repetition starting from the final state of the previous run. A small number of 
such repetitions is sufficient to reduce the scaling of the mixing time M(e) to 
logarithmic in e. 

Both the quantum walk on the line and the cycle thus provide a quadratic 
speed up over classical random walks. This quadratic speed up does not carry 
over to all quantum walks on higher dimensional structures, see, for example, 
IE M, M ■ It remains an open question how ubiquitous this behaviour really 



is. 
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3 Non-unitary quantum walks 



The quantum walk dynamics described in the previous section is a pure quan- 
tum evolution terminated by a final measurement to determine the outcome 
of the quantum walk. If, instead, the quantum walk is measured after every 
step, it is easy to see that it becomes a classical random walk. In between 
these two extremes, a smaller number of measurements can be included in the 
quantum walk dynamics, as in the "warm start" already mentioned for the 
walk on the cycle. This is not a new idea. Quantum walks with measurements 



were first explored in |2l| for controlling a physical quantum walk in an atom 
optical system. 

We are thus led to define a more general quantum walk by including the pos- 
sibility of non-unitary operations at each step (pure quantum dynamics is uni- 



tary, measurements are non- unitary). Kendon and Sanders [22| provide a gen- 
eral formulation of non-unitary quantum walks using superoperators. We will 
not need the full generality here, we restrict this study to randomly-occurring 
uncorrelated non-unitary events, and measurements at specific chosen times 
during the quantum walk. Nonetheless, the evolution of the quantum walk 
must now be described using a density operator p(t) given by 

p(t + l) = (l-^SCp^C^+p^TPj-SCp^C+StpJ. (12) 

3 

Here Pj is a projection that represents the action of the non-unitary operator 
and p is the probability of applying this operator per time step, or, completely 
equivalent mathematically, to a weak coupling between the quantum walk 
system and a Markovian environment with coupling strength p. For a pure 
state, the density operator p(t) = |\E'(t))(\I/(t)|, it thus has the normalisation 
Tr[p(t)] = l. 

Full analytical solution of a non-unitary quantum walk has been done only 
for a few instances. For quantum walks on the line, Brun et al analysed 



the case of random measurements on the coin only. While analytical solution 
is challenging, eq. (12) lends itself readily to numerical simulation since p, 
S and C can be manipulated as complex matrices, while the Fj generally 
remove some or all of the off-diagonal entries in p. Kendon and Tregenna 



121 ] evolved eq. (12) numerically for various choices of Pj, projection onto the 
position space, projection into the coin space in the preferred basis (| ± 1)), 
and projection of both coin and position. In all cases, the spreading rate is 
reduced, in the long time limit 23], it becomes proportional to yT instead of 
proportional to T. More interesting behaviour is seen for intermediate times 
and noise rates p with noise applied to the position. Kendon and Tregenna 
observed that, for 2 < pT < 5, the distribution becomes very close to uniform 



while retaining the full quantum linear spreading rate [12] . With noise applied 
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to the coin only, the distribution retains a cusp shape. To quantify this the 
TVD given by eq. (9) is used, this time with P u defined to be a top-hat of 
appropriate width x £ {±T/y/2}, see [s], 0]. 

For a quantum walk on the cycle subjected to Markovian noise, the mixing 
behaviour is dramatically improved, provided noise is applied to the position. 
The noise guarantees mixing to the uniform distribution, and a similar judi- 



cious choice of 2 < pN < 5 produces the minimum mixing time [12J. Noise on 
the coin only does cause the quantum walk on a cycle to mix, but not signifi- 
cantly faster than a classical random walk. Furthermore, the time-averaging is 
no longer needed, fast mixing occurs in time 0(iV log(l/e), as shown in numer- 



ical work by Maloyer and Kendon [15|, [24J , and recently proved by Richter [14 



Significantly, Richter also proved that the randomness produced by quantum 
measurements alone is sufficient to produce an optimal mixing time. By apply- 
ing measurements at regular intervals, instead of randomly with probability 
p, the speed up is still obtained. 

We thus have two examples where the optimal computational properties are 
obtained for a judicious combination of quantum dynamics and measurements, 
which introduce a component of randomness into the otherwise deterministic 
quantum dynamics. In order to understand more about how this mixture 
of quantum and random resources combine to produce their computational 
power, we next describe how we can compare them in the quantum walk. 



4 Comparing quantum and classical computational resources 



Our goal is to compare the resources required to perform a quantum walk 
with those required for a classical random walk, to find out which is more 
efficient when the random resources are taken into account as well as the 
memory and gate operations. What we will find is that the optimal non- 
unitary quantum walk on the line run for T steps requires only O(logT) extra 
quantum gates and ancillae compared with a pure quantum walk. Since the 
pure quantum walk uses O(TlogT) quantum gates and O(logT) qubits, the 
extra resource costs are insignificant to leading order. A classical random walk 
producing the same computational outcome requires T 2 steps. This can be 
accomplished using 0(T 2 logT) quantum gates and 0(T 2 ) ancillae. The extra 
resources compared with a fully quantum walk are mainly needed to generate 
the random bits. Applying the same reasoning to the optimal quantum walk 
on a cycle shows that it is similarly efficient. We now explain in detail how 
these results are obtained. 

We can quantify the randomness added to, or generated by, the non-unitary 
quantum walk by counting up the number of random bits, either supplied 
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to time the random measurements, or produced from the results of those 
measurements and then used as input to the next iteration of the quantum 
walk. The quantum resources can be quantified by counting the number of gate 
operations required, in a similar manner to classical computational complexity. 
What we then need is a way to compare the quantum and random resources. 

We first observe that, if we have only quantum resources to hand, we can 
generate random bits by using a quantum gate followed by a measurement. 
To illustrate, using a quantum coin by itself, we apply the coin toss then 
measure. For example, 

C| + l) = {|-l) + |+l)}/v / 2. (13) 

The measurement outcome is | — 1) 50% of the time, and |+1) the other 50% 
of the time, i.e. one random bit. We can repeat this procedure (starting with 
the coin in whichever state was obtained after measuring) to obtain further 
random bits. Of course, the converse does not work: if we have only classi- 
cal random bits, we cannot efficiently simulate all the properties of a general 
quantum system. What is notable is that a quantum computer comes with 
randomness as a "built in" function, whereas classical computation does not 
and it must be added separately, or faked (pseudo-randomness) at some sig- 
nificant computational cost. 

An alternative way of viewing non-unitary quantum evolution is as a pure 
quantum evolution in a larger system that includes an ancilla or environment 
degrees of freedom. Let us illustrate this with the quantum coin in an arbitrary 
state 

|# c )=a|-l)-/?|+l). (14) 

If we measure in the computational basis, we find the coin in state |— 1) with 
probability a 2 and in state | + 1) with probability f3 2 . Now we add an ancilla 
qubit, another two state system with states |eo) and \e\). We start with the 
ancilla in state |e ) and apply a CNOT gate with the quantum coin as the 
control and the ancilla as the target, 

CNOT|* c )|e ) =a|-l)|e ) -/3|+l)|e!). (15) 

Now we discard the ancilla, and consider the state of the quantum coin only. 
Mathematically, we trace out the ancilla, leaving 

Tt e [a|-l)|e > - /3|+l)| ei >] = q; 2 |-1><-1| + /3 2 |+l)(+l|. (16) 

This is now a mixed state with classical probabilities that match the proba- 
bilities measured in the quantum state eq. (14). However, it now has no useful 
quantum properties. Notice in particular that the phase {—(3 rather than +/3) 
is lost from the original superposition state. The quantum phases are cru- 
cial for generating the computational speed up in a quantum walk, via the 
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cancellations (interference) they produce when the quantum walk arrives at 



the same location via two different paths. Kendon and Sanders |22| employ 
this method of measuring the quantum coin with a more general coupling to 
explore the effects of weak measurements. 

Based on these examples, we can see that one way to count resources in a non- 
unitary quantum walk is to count them all as quantum resources. We count 
the number of quantum gates, and we also tally up the number of quantum 
bits (qubits) we require, including the ancillae that we use in the measurement 
process. A random bit thus requires one quantum gate, and one qubit ancilla 
to couple to the qubit being measured. Note that we cannot "recycle" the 
ancillae, since returning them to a pure quantum state \bq) would require 
further operations and ancillae. 

We next apply the quantum resource counting to the quantum walk on the line, 
run for T steps with non-unitary noise rate p. This is equivalent to including 
a total of pT measurements at random intervals during our quantum walk. 
Since the noise can be applied as a weak coupling at every step, we do not 
need to generate random numbers to decide when to apply the measurements. 
Each step of the quantum walk uses two quantum gates, C and S. Now, C 
acts non-trivially on the coin qubit only, but S acts on the whole quantum 
walk system. We assume the position of the quantum walker is encoded in a 
quantum register of size |Tog(2T + 1)] qubits. Thus S will need to consist of 
O(logT) elementary gates each acting on one or two qubits. Next we turn to 
the noise. There are two cases we need to consider: noise acting on the coin 
only, and noise acting on the position space as well. (Noise acting on both the 
position and coin differs only 0(1) in resources required from noise acting on 
the position only, so we will consider these two cases as one.) We have already 
noted that noise on the coin only does not have the same computational effect 
as noise acting on the position. We will now see that the required resources 
are different. The pT random measurements applied to the coin only require 
0(pT) quantum gates, and 0{pT) ancillae. Applying noise to the position 
requires 0(pT log T) quantum gates and 0(pT log T) ancillae. These random 
bits are, of course, drawn from the distribution of the quantum walk on the 
line, which is not uniformly random, so this estimate is an upper bound on 
the true minimum required 1 1. The number of qubits we need for the quantum 
walk is one for the coin plus [log(2T + 1)] for the position. We thus require a 
total of 

coin noise position noise 

quantum gates: O (T log T + pT) O (T log T + pT log T) 

qubits: 0(logT + pT) 0(logT + pT logT) 



1 The distribution of the quantum walk on the line is fairly close to uniform [if, 
so this will be a good estimate of the requirements 
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Let us check the extreme cases, pure quantum and fully classical. For the 
pure quantum walk (no noise or randomness) we have p = 0, so for large T 
we require O(TlogT) quantum gates and O(logT) qubits. For the classical 
random walk, p = 1 corresponding to measurements applied at every step, but 
it is sufficient to apply measurements to the coin only, so we find we still require 
O(TlogT) quantum gates, but we need an additional T ancillae to generate 
the randomness at each step. However, to make a fair comparison in terms 
of the outcome of the walks, the classical random walk must run for longer. 
To achieve the same spreading as a quantum walk run for T steps, we require 
~ T 2 /2 steps of a classical random walk, requiring 0(T 2 logT) quantum gates 
and 0(T 2 ) ancillae. Viewed from this perspective, randomness is an expensive 
resource! Of course, what we did not count is the cost of maintaining the 
quantum coherences necessary for the pure quantum walk to function properly. 
However, this depends on the particular physical implementation of a quantum 
computer. While there is no fundamental minimum requirement that we know 
of, in practice we expect the costs of maintaining coherence in a realistic 
quantum computer to be very high in terms of the quantum error correction 
overhead that will be required [251 ] . 

Finally we are ready to consider the optimal quantum walk. As already noted, 
the closest to uniform distribution requires measurements to be applied to the 
position, and is obtained for a noise rate such that pT ~ 0(1). Our resource 
requirements under this condition are 0(TlogT + logT) quantum gates and 
O ((1 + 0(1)) logT) qubits. We thus need to add only O(logT) extra quantum 
gates and ancillae, to optimise the quantum walk, a small increase in the 
resources compared to the total resources required for the pure quantum walk. 

We can do the same comparisons for the non-unitary quantum walk on the 
cycle. To obtain a distribution within e of the uniform distribution on a cycle 
of size N, we run for M(e) steps with noise rate rate p. 

coin noise position noise 

quantum gates: 0(M(e) logiV +pM(e)) 0(M(e) logiV + pM(e) logiV) 
qubits: 0(log N + pM(e)) 0(\og iV + pM(e) log N) 



For the time-averaged mixing time M(e), we need to add 0(logM(e)) gates 
and qubits, to choose a random time at which to stop between . . . T, as per 
the definitions in eqs. (10) and (11). 

Using M(e) ~ 0(iV log(l/e)) for the quantum walk with noise rate p such 
that pN ~ 0(1), and M(e) ~ 0(N 2 log(l/e)) for a classical random walk, we 
can compare the resources required in each case. Noise on the coin applied at 
every step (p = 1) is sufficient to produce a classical random walk, which thus 
requires 0(N 2 log iVlog(l/e)) quantum gates and 0(N 2 log(l/e)) qubits. For 
the optimal quantum walk we consider only noise on the position, since noise 
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on the coin only does not allow the quantum walk to mix significantly faster 
than a classical random walk [l2j]. A pure quantum walk run for iVlog(l/e) 
steps (which does not mix) requires 0(iVlog iVlog(l/e)) quantum gates and 
O(logiV) qubits. For optimal mixing we need to add 0(logiVlog(l/e)) quan- 
tum gates and qubits. Again, only a small amount of randomness is required to 
optimise the quantum walk on a cycle, while the cost for the classical random 
walk is much higher. 



5 Entanglement in mixed states 



Having estimated the number of quantum gates required for the quantum 
walk, it would be interesting to know whether they are actually used to full 
effect. The unitary operations that generate quantum correlations are not 
guaranteed to do so every time they are applied. Depending on the current 
state of the system, they can just as easily remove correlations as add them. 
So, we will look directly at the quantum correlations in the quantum walk 
system by calculating the entanglement between the coin and the position of 
the quantum walker. We choose our entanglement measure to be the nega- 
tivity because this can be calculated numerically in a fairly straightforward 
manner for density operators such as p{t), and there are few options that 
meet this criterion. This is not ideal. The negativity does not relate directly 
to information-theoretic quantities, like entropy, that could be compared with 
the number of random bits, but it will suffice for our purposes in this study. 

The negativity is defined as follows. First we must choose a division of our 
system into two (or more) subsystems between which to identify the entan- 
glement. For our quantum walk, the natural division is between the coin and 
the quantum walker's position. We note that the entanglement across this 
division will be the same whether we regard the quantum walk as a qubit 
coin and a unary position, or as the position encoded in a binary quantum 
register (because the Hilbert spaces are the same size). We perform a partial 
transpose on one subsystem to obtain a new matrix p'(t). For example, the 
partial transpose with respect to the coin subsystem is 

Pxcybtt) = Pxb,yc(t) (17) 

where x, y are position indices and c, b are coin state indices. Next, we de- 
termine the spectrum of p'(t), denoted by{A^}. The normalisation of pit) is 
carried over to p'(t), so — 1; but unlike pit), it is possible for p'(t) to 



have negative eigenvalues. The negativity is defined [26|, |27|, [28| as 

E 



H? |A H 



10 



time step 



Fig. 1. Negativity for noise rate p = (circles), p = 0.05 (asterisks) and p = 0.1 
(squares) for a quantum walk on the line with noise applied to both coin and 
position, run for 50 time steps. 

which is just the sum of the negative eigenvalues. The negativity ranges 
between zero and one, with any non-zero value indicating entanglement is 
present. If the negativity is zero it means the state is probably not entangled, 
but there can be exceptions 29[. The exceptions are known to be relatively 



rare in the set of all possible states [30|, and the entanglement they contain 
is difficult to apply to useful quantum tasks [2^]. For this study, we will not 
need the fine-grained detail of these possible exceptions. 



6 Entanglement in a quantum walk on the line 



The entanglement in a pure state quantum walk has been studied previously, 



see for example, [31]. It fluctuates with each step, and eventually settles down 
to an asymptotic value that depends on the initial state of the quantum coin, 
and on any bias in the quantum coin operator C. The addition of noise smooths 
out this behaviour, see fig. 1, and steadily reduces the level of entanglement 
between the coin and the position. Recall that a noise rate p can be inter- 
preted as a weak measurement of strength p applied every step, or, a full 
measurement applied only with probability p, whereupon the steady reduc- 
tion in entanglement can be regarded as the average effect of a few random 
measurements. Our simulations simply take eq. (12) and evolve it numerically, 
calculating the negativity and TVD using an appropriate top hat distribution, 
for various types of noise and noise rate p. 



Guided by the results in [12j and [31[ , we focused on the entanglement at the 
end of the quantum walk, just before the final measurement, and considered 
how it varies as the noise rate p is varied. We examined three cases of noise, 
applying it to both the coin and position, and also separately to just the coin 
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Fig. 2. Negativity (diamonds) and TVD for noise applied to position (dashed), coin 
(dotted) and both (solid) for a quantum walk on the line of 100 steps. 

or the position. An example is shown in fig. 2. Both the TVD from the optimal 
top hat distribution, and the negativity are plotted. The minimum in the TVD 
indicates the optimal measurement rate p. For a quantum walk of 100 steps 
the optimal p ranges through 0.025 < p < 0.05. Note that the minimum for 
measurements on the coin only is much shallower. The distribution in this 
case retains a cusp shape [12|. The negativity also remains above zero for 
longer than when noise is applied to the position of the walker, indicating the 
coin measurements are less effective at removing the quantum correlations. 
If noise is applied to the position, with or without noise on the coin as well, 
the negativity drops to zero at p ~ 0.055, shortly after the optimal noise 
rate is reached. The optimal amount of randomness is thus just about the 
amount required to remove all the the quantum correlations from the system. 
This make intuitive sense: we are trying to achieve a uniform distribution 
on the line, in which any location within the top hat region is equally likely. 
Quantum correlations distort this smooth distribution, giving it peaks and 
troughs, especially at the ends of the top hat 0, @]. If the noise rate is turned 
up until the classical random walk is obtained for p = 1, classical correlations 
build up to produce the binomial distribution in which the quantum walker is 
more likely to be found nearer the starting point of the walk. 



7 Entanglement in a quantum walk on cycles 



For the mixing time on cycles, there are extra considerations. The mixing time 
in general depends not only on the size of the cycle, but also on how close 
to uniform one sets the threshold e. As already noted, pure quantum walks 
don't mix unless something is done to disrupt the pure quantum evolution, 
and random or even regular repeated measurements efficiently change the 
behaviour into that of fast mixing to the uniform distribution. Note also that 
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time step 



Fig. 3. Negativity (diamonds) and TVD for noise applied to position for a quantum 
walk on the cycle of size 29 with p = 0.2239 (dotted), p = 0.2511 (solid) and 
p = 0.2818 (dashed). The grey horizontal line is at 1/29. 



the optimal rate of measurement 2 < pN < 5 is independent of the threshold 
e, even though the main effect is to provide logarithmic scaling of the mixing 
time with e. 



We studied how the entanglement varies during the quantum walk on a cycle, 
with the noise rate p chosen to be near-optimal. As with the walk on the line, 
we examined the three cases of noise applied to both the position and coin, and 
applied to just the position or coin separately. Again noise applied to the coin 
only does not provide a significant improvement in the behaviour compared 
with noise applied to the position. The results for a typical example, a cycle of 
size A = 29 with noise applied to the position only, are shown in fig. 3. While 
the actual mixing time is determined by the choice of e, we have indicated 
the position of 1/N in fig. 3, this being the probability of finding the walker 
at one location in a uniform distribution. The time at which the TVD drops 
below this line is around the time the entanglement also drops to zero. The 
variability of the TVD shows that the mixing time is not a smooth function 
of e, so we cannot expect to determine a more precise result. As we argued for 
the walk on the line, quantum correlations are incompatible with the result 
we want, a uniform distribution, so the optimal quantum walk arranges for 
the quantum correlations to be removed by the end of the process. 



13 



8 Discussion 



The optimal quantum walk on the line and cycle is a carefully balanced combi- 
nation of quantum dynamics with randomness provided by repeated measure- 
ments during the evolution of the walk. We have shown that, while the quan- 
tum correlations are necessary to obtain linear spreading and mixing times, 
they must be neutralised to produce a uniform final distribution. We have also 
analysed the resources required for a non-unitary quantum walk and shown 
that the additional randomness required to obtain the optimal behaviour is 
only logarithmic in the size of the problem. Since a quantum computer has 
randomness as a built in function, this is a very efficient use of resources com- 
pared to performing a classical random walk, which requires an amount of 
randomness quadratic in the size of the problem. 

While the relationship between quantum entanglement and shared random- 



ness is well-studied in quantum communications theory, for example, see [32 
the relationship between the two types of resources in quantum computa- 
tion has received little attention. The idea of combining measurements with 
quantum gates is well-established [33], but little is known about the impor- 
tance of each contribution to the power of the computational process. Note in 
particular that in the cluster state model, random measurement outcomes are 
corrected, requiring extra computational steps, whereas the non-unitary quan- 
tum walk takes advantage of this randomness to enhance the computational 
efficiency. Our results contribute to the task of quantifying this relationship 
and to furthering our understanding of how to harness the computational 
power of quantum systems. 

We thank many people for interesting discussions of quantum walks, among 
them, Ivens CarneircT 2 ! Hilary Carteret, Jochen Endrejat, Barbara Kraus, 
Peter Richter, Barry Sanders, Mario Szegedy, and Tino Tamon stimulated 
our thinking for the work in this paper. VK is funded by a Royal Society 
University Research Fellowship. 
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